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ABSTRACT 

A number of authors have argued that the Sun must have been born in a cluster of no more than 
several thousand stars, on the basis that, in a larger cluster, close encounters between the Sun and other 
stars would have truncated the outer Solar System or excited the outer planets into eccentric orbits. 
However, this dynamical limit is in tension with meteoritic evidence that the Solar System was exposed 
to a nearby supernova during or shortly after its formation; a several thousand star cluster is much 
too small to produce a massive star whose lifetime is short enough to have provided the enrichment. 
In this paper we revisit the dynamical limit in the light of improved observations of the properties 
of young clusters. We use a series of scattering simulations to measure the velocity-dependent cross- 
section for disruption of the outer Solar System by stellar encounters, and use this cross-section to 
compute the probability of a disruptive encounter as a function of birth cluster properties. We find 
that, contrary to prior work, the probability of disruption is small regardless of the cluster mass, and 
that it actually decreases rather than increases with cluster mass. Our results differ from prior work for 
three main reasons: (1) unlike in most previous work, we compute a velocity-dependent cross section 
and properly integrate over the cluster mass-dependent velocity distribution of incoming stars; (2) we 
recognize that ~ 90% of clusters have lifetimes of a few crossing times, rather than the 10 — 100 Myr 
adopted in many earlier models; and (3) following recent observations, we adopt a mass-independent 
surface density for embedded clusters, rather than a mass-independent radius as assumed many earlier 
papers. Our results remove the tension between the dynamical limit and the meteoritic evidence, and 
suggest that the Sun was born in a massive cluster. A corollary to this result is that close encounters 
in the Sun's birth cluster are highly unlikely to truncate the Kuiper Belt unless the Sun was born in 
one of the unusual clusters that survived for tens of Myr. However, we find that encounters could 
plausibly produce highly eccentric KBOs such as Sedna. 

Subject headings: planets and satellites: formation — solar system: formation — stellar dynamics 



1. INTRODUCTION 

It is very likely that our Sun, like most stars, was born 
as part of a star cluster. The birth cluster environment is 
a self-gravitating system of gas and dust that collapses at 
a multitude of points to create its membership of stars. 
This process is inefficient, converting only 10% to 3 0% of 



the mass to st ars while leaving 70% to 90% as gas ( |Lada 



& Lada 2003). As long as the gas remains, the stars are 
gravitationally bound and orbit chaotically throughout 
the cluster. This motion continues until the winds and 
ionizing radiation of the newborn stars blow away the 
remaining gas. At that point the cluster is no longer 
self-gravitating, it disperses, and the stars are free to 
move about the Galaxy. 

By examining the state of our Solar System as it exists 
today, we can place limits on the birth cluster's prop- 
erties (i.e. mass, surface density, radius, etc.) and thus 
have a better understanding of the environment in which 
our Solar System was born. 

Meteorites provide one line of evidence. Samples taken 
from meteorites that solidified not long after the Solar 
System formed show daughter products of multiple ra- 
dioactive isotopes. These radioactive isotopes most likely 
originated from stellar nucleosynthesis, and their pres- 
ence suggests that there was a nearby supernova that 
enriched the Solar Nebula dThrane et al.|2006j |Wadhwa| 
20071 lAllen et al.l|2007l [Connelly et al.||2U08| |Duprat &| 
Tatischefi||2008| |Adams||2Q10p . Such early exposure is 



likely only if the birth cluster includes a star > 25M in 
mass. This consider ation suggests a birth cluster with at 
least 1000 members ( Adams|2010 ). A priori such a large 
cluster is not unlikely. The young cluster mass function 
is 



(IN oc Af~ 2 over the range from ~ 10 1 - 1Q 6 M K 

hich 



Lada 



|fc Lada|[2003l |Chandar||200^ |Fall et aL|[2010] ) , wine 
plies that ~ ^ — | of stars are born in clusters larger than 
this0 

On the other hand, several authors have attempted to 
obtain upper limits on the birth cluster size by consid- 
ering the effects of close encounters with other stars on 
the young Solar System. Close encounters would perturb 
the Solar Nebula during its formation. Any such pertur- 
bation must still permit the formation of a Solar System 
like the one we see today, with eight planets that all lie 
close to the same plane (inclination angles < 3.5°) and 
have small eccentricities (less than 0.2, or less than 0.09 

1 We do caution here that the observations do not completely 
establish the existence of a mass function dN/alM oc M~ 2 over 
this full mass range. Th i s pow er law is seen in both the Galactic 
sample o f jLada fe L ada (2003| and the extragalac tic samples of 
iChand ar (2UU9) and Fall et al. (2010}. However, the jLada fe Lada| 
( |2UU3| ) sample only includes clusters wit hin 2 kpc ot the bun, the 
most massive of which is ~ 10 3 Mq. The |Chandar| ( |2QQ9] ) and |Fall| 
|et al.| ( |2010| samples only go down to a mass of ~ 10 4 Mq . Thus 
the data do not rule out the possibility that there is a break in the 
mass function in the range 10 3 - 10 4 Mq. However, it would be 
quite a coincidence for this discontinuity to coincide exactly with 
the mass range where the observations are incomplete. A simpler 
hypothesis is that dN/alM oc M~ 2 over the full mass range. 
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is we exclude Mercury, whose eccentricity is pumped by 
p erturbations from other p lanets) . 



Adams & Laughlin ( 2QQ1[ ) determined an upper bound 
on the birth cluster size by performing Monte Carlo scat- 
tering simulations of close encounters between the Solar 
System and passing stars. They found that the maximum 
number of stars the birth cluster could have had without 
causing a Jovian planet to be excited to an eccentricity 
greater than 0.1 is 2200 ±1100. However, observations of 
the properties of young, embedded star clusters were at 
the time quite limited, and they were forced to make a va- 
riety of assumptions about cluster properties. First, the 
relative velocities of the stars in their simulations were 
chosen randomly from a Maxwellian distribution with a 
standard deviation of 1 km/s; thus they had very few 
simulations done with higher velocities, and the velocity 
distribution was implicitly assumed to be independent of 
the cluster mass. Second, they assumed that the clus- 
ter's lifetime, and thus the time the Sun was exposed to 
close encounters, was on the order of 100 cluster crossing 
times. Lastly, they assumed that the cluster's radius was 
fixed and independent of its mass. This last assumption 
would lead a high mass cluster to have very high surface 
densities. 

More recent observations have shown that most of 
these assumptions are not satisfied in typical embedded 
clust ers. Compiling obser v ations of embedded cl usters 
from |Shirley et al. (2003), Faundez et ah] fl2004p , and 
Fontani et al | ( |2005| ) ^Pall et al.| ( |2010j ) hnd thaTembed- 
ded clusters form a sequence ol roughly constant surface 
density, not constant radius. Similarly, current observa- 
tions f avor cluster lifetimes closer to one to four crossing 
times (|Elmegreen||2000[ |Hartmann et al.] 2001 |Lada fc 
Lada|2003| I'l'an et al.||2006|| Jeffries et al.||2011[ |Keggiani 
et al ]|201ip , not 100. Mter ~ 10 Myr, 9QWot the stars 
born m clusters have dispersed into the field, and by 
100 Myr roughly 99% have done so fl Allen et al.||2007 



Fall et al.||2009| Chandar et al.||2010[ ). We do not know 
the lifetime ol the Sun's birth cluster, but in obtaining 
limits on its size it is important to consider the possibil- 
ity both that the lifetime was very short, and that it was 
tens of Myr, rather than assuming the latter as was done 
in most previous work. 

In this paper we show that retaining the sa me theoret- 
ical framework as Adams & Laughlin (2001), but using 
these more recent observations to set the properties of 
the embedded clusters leads to a dramatic relaxation of 



the Adams & Laughlin (2001) limit on the size of the 
Sun s birth cluster. In turn, this greatly eases the ten- 
sion between the meteoritic and dynamical constraints. 
In particular, the lifetime of a 25 M star before it goes 
supernova is ~ 7.54 Myr, longer than most clusters sur- 
vive. A 75 Me?) sta r, however, has a lifetime of ~ 3.64 
Myr. Adams (20101) finds that to have a reasonable prob- 



rent picture of Solar System formation. The existence of 
Neptune, with its size and orbit, suggests that the surface 
density of the Solar Nebula should have been high enough 
to produce Kuiper Belt Objects (KBOs) out to distances 
greater than 50 AU with a gradually declining object den- 
sity. An unperturbed Solar Nebula, which tends to pro- 
duce objects with low eccentricities, is hard to reconcile 
with Sedna's eccentricity of 0.84 and semi-major axis of 
542 AU. A close encou nter with another star could pro- 
duce a disk truncation ( Ida et aLp OOO ; Ko bayashi fc Ida| 
2001) and scatt er objects into high eccentricity orbits to 
produce Sedna (Brasser et al.||2006| ). 

Close encounters could also affect other putative ob- 
jects orbiting th e Sun: the proposed companion brown 



dwarf Nemesis rtRaup fc Sepkoski Jr. 1984 Whit mire 
fc Jackson IV 1984) or the more recently proposed Ty- 



che"| Matese 1999 nMatese et al.||2005l |Matese fc Whit 
|mire||201l[ ). Tyche has been proposed to have an orbit 
just inside the Inner Oort Clound (a semi-major axis of 
5,000 AU) with a mass of 1-4 Juipter Masses. G iven the 



low p robability of capturing such a companion (Tohline 
2002[ ), if either Nemesis or Tyche exists, they must have 



formed with the Sun. To date, there has been no re- 
search on whether Nemesis nor Tyche could survive the 
dynamics of the birth cluster. 

In this paper we address all of these issues. We have 
con struct ed simulations similar those of |Adams fc Laugl>] 
lin (2001). We use these together with updated obser- 
vational data on the properties of embedded clusters to 
determine an upper limit of our birth cluster's member- 
ship size and the effects of a close encounter on both 
the truncation on the Kuiper Belt and production of a 
Sedna-like object. In addition, we will use this data to 
determine the likelihood of Tyche staying bound to the 



Sun d uring its time within the birth cluster. | Adams et al. 
(2006) have undertaken similar calculations, but survey 



ability ol producing such a star, a birth cluster of mem- 
bership size > 10 4 is needed. Again, such a large cluster 
is not a priori unlikely, given the observed young cluster 
mass function. 

Close encounters within the Sun's birth cluster have 
also been theorized to be responsible for other aspects 
of the Solar System's present-day architecture, particu- 
larly the drastic drop in object density in the Kuiper Belt 
at ~45AU, and the Kuiper Belt object Sedna. Both of 
these phenomena are difficult to understand in our cur- stellar cluster is Maxwellian (Proszkow & Adams 2009) 



a much more limited part of parameter space. In par- 
ticular, they do not consider clusters of more than 10 3 
stars. In Section 2 we will look at the simulation and 
method used in our experiment. In Section 3 we will 
then analyze the results of the simulation and finally in 
Section 4 we will make our conclusions. 

2. METHODS 

To determine an upper limit on the birth cluster mem- 
bership size, we want to determine the probability of a 
disruptive encounter as a function of cluster surface den- 
sity and mass. By disruptive event, we mean a close 
encounter in which any of the Jovian planets are excited 
to an eccentricity greater than 0.1. In the process, we 
will also determine whether encounters with passing stars 
could be responsible for truncating the Kuiper Belt at 
~50AU, for producing a Sedna-like object, or for strip- 
ping any distant brown dwarf companion. 

2.1. Determining Probability 

To obtain the probability of an event occurring, we 
first perform numerical scattering simulations to deter- 
mine the velocity-dependent cross-section for that event, 
(a). We then calculate the average cross section times ve- 
locity, (crv). We defer a description of our simulations to 
section 2.2. Simulations of star clusters suggest that the 
distribution of encounter velocities within a embedded 



3 



so 



(av) 



(a)v 3 e 2ct - dv 



(i) 



where a v is the standard deviation of the velocity distri- 
bution. If we then multiply this by the stellar number 
density of the cluster, n, we obtain the encounter rate, 



T = n(av). 



(2) 



If we multiply T by the total time the Sun was exposed in 
the cluster, we will have the expected number of events, 
A. This time will be roughly the time for which the clus- 
ter survives before mass evaporation leads it to disperse, 
tnf e . Thus, we have 



A = n life . 



(3) 



To compute the overall probability that an event will 
occur, we will assume that events follow a Poisson distri- 
bution, in which case the probability that an event will 
occur at least once is 

Pevent ~ 1 ^ • (4) 

The calculation is slightly more complicated for the 
Kuiper Belt, where there are many KBOs, and vary- 
ing numbers can be affected by a single encounter. In 
our simulations we represent particular radial bins of the 
Kuiper Belt with a number N t of test particles. We use 
our simulation to determine (a)k, the velocity-dependent 
cross section for an event (e.g. excitation or ejection from 
the Solar System) to affect k of these N t test KBOs. Us- 
ing equation]!] we can then calculate (crv)k for each value 
of h. The rate at which the KBO population undergoes 
events is then given by 



N t k 

k=i 1 



(5) 



where n is the number density determined in equation 
[7[ Once V has been determined, we use equation p| to 
compute A. The expected fraction of KBOs that undergo 



compute A. The expected 
a certain even is then 

</> 



l-e- A . (6) 

The initial conditions of the birth cluster are unknown 
to us, so we cannot directly compute all of these quan- 
tities. However, we can estimate them in terms of the 
cluster surface density, £ c , and cluster mass, M c . For a 
spherical cluster, the number density, velocity dispersion, 
and crossing time are 



3^/2^3/2 



4raM, 



1/2 



700M~ 1/2 ^ pc~ 3 



1/2 



= 3.2^( r 2 Ml /4 E 1 _ / 1 4 km s- 1 
Ml" 



(7) 



(8) 



G 1 / 2 (7rS c ) 3 / 4 



0.62M 4 1/4 S_i /4 Myr, (9) 



where a v i r is the virial parameter, m is the mean s tel- 
lar mass, which we take to be O.2M ( Chabrier|2QQ5| ) for 



our calculations, M4 = 
2 Equations (|7|) 



M c /10 4 M , and E_i = E c /0.1 
g cm - ". Equations (IYJ) and (ISJ) provide estimates of the 
number density and velocity dispersion that we can use in 
our estimates. The relationship between the cluster life- 
time tiif e and the crossing time, however, is somewhat un- 
certain. As discussed in the Introduction, most clusters 
disperse rapidly, although exactly how rapidly is debated. 
For example, the best studied cluster containing O stars 
is the Orion Nebula Cluster [M 4 = 0.46, E_i = 1.2 - 
Hillenbrand fc Hart mann|| 1998] ), for which t cross = 0.44 
Myr. Es timates ol its "age spread range from ~ 1 — 3 



Myr. Estimates 01 its a ge spread range Iro m ^ 1 — d 
Myr (e.g. [Reggiani et al.|[20Tl| | Jeffries et al ||2011[ ), cor- 



responding to ~ 1 — 6£ croS s 5 depending on the technique 
used to estimate protostellar ages and similar systematic 
uncertainties. On the other hand, ~ 10% of clusters sur- 
vive for more than 10 Myr. We therefore consider two 
possible scenarios for the lifetime. The more likely one 
is that the Sun was born in one of the rapidly-dispersing 
clusters, and in this case we take t\a e = t cross ; the true 
lifetime may be a factor of a few larger, but, as we will 
see below, this factor of a few is relatively insignificant. 
The other possibility is that the Sun was born in the 
minority of clusters with a longer lifetime. To represent 
this case, we take t\a e = 30 Myr, typical of this class. 

Finally, note that, in the rapidly-dispersing cluster sce- 
nario, for massive and high surface density clusters the 
cluster lifetime may be significantly less that the main 
sequence lifetime even of a very massive star. This is 
certainly the case in the ONC, for example, where star 
formation is complete and the cluster is likely in the pro- 
cess of dispersing, but 9 1 Ori C is still alive. This is not 
necessarily a problem for enrichment, however. Even for 
a dispersing cluster, the relatively low velocity disper- 
sion of the stars implies that stars will not spread that 
far before the most massive stars go supernova. 

2.2. Simulations 

The simulations we use to measure (a) consist of the 
Solar System and the scattering star system. The So- 
lar System includes the Sun, the four Jovian planets, 
Tyche and 320 KBOs all within the same plane. We 
only place Tyche in the simulation because the semi- 
major axis of Nemesis is so large that we cannot rea- 
sonably evaluate it behavior with scattering simulations. 
Instead, a full (and much more expensive) N-body simu- 
lation of the entire birth cluster would be required. We 
treat the Jovian planets and KBOs as massless to con- 
serve computing power. This means that we are unable 
to model secular processes produced by planet-planet in- 
teractions. However, since such processes generally take 
much longer than that ~ 0.01 Myr covered by our sim- 
ulations, we could not model these in any event. We 
consider Tyche masses of zero, 0.02 M and 0.08 M . 
We set the orbital properties of the Jovian planets to 
their present-day values and we randomize the initial or- 
bital phases. We give Tyche a semi-major axis of 1,000 
AU^J an eccentricity randomly between 0.0 and 0.4, and 
a random argument of periastron. We set up the Kuiper 

While t his i s closer than the 5000 AU proposed by |Matese| 



|et al.| ( pb05| > and |Matese W hitmire (2011), the use of a similar 



but smaller distance greatly reduces the computational cost by 
allowing us to consider a smaller range of incoming star impact 
parameters. 
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Fig. 1. — The log of the probability of a close encounter exciting each Jovian planet's eccentricity to greater than 0.1, in one crossing 
time, as a function of cluster mass M c and surface density E c . 



Belt up such that all objects are placed in 32 concentric 
circles ranging from 35 AU to 500 AU, in 15 AU inter- 
vals. Each circle has ten objects equally spaced in angle, 
and all objects have a zero eccentricity. For our statisti- 
cal analysis, we then divide the KBOs into eight distance 
bins (35-80AU, 95-140AU, 155-200AU, 215-260AU, 275- 
320 AU, 335-380 AU, 395-440AU and 455-500AU), each 
bin containing four concentric circles. 

The incoming star systems consist of either a solo star 
or a binary star system. We first select the mass of 



drawing if the mass exceeds the primary's mass. Finally, 
we select the orbital period of the binary by drawing 
from the distribution obser ved for field G star binaries 
flDuquennoy fc Mayor]|199lD 



p(logr) 



1 



2W 



(11) 



the p rimary incoming star by drawing from a Chabrier 
( |2005| ) IMF. We then determine if the primary has a 
companion. We set the probability of there being a com- 
panion to 



fi 



binary 




: m < 0.5 M© 
: 0.5 M < m < 2.0 M 
: m > 2.0 M 
(10) 

based on a rough fi t to the mas s-dependent companion 
fraction reported by |Lada (2006). If a companion exists, 
we also select its mass from the IMF, rejecting and re- 



where r is the period in days, a\ ogT =2.3, and logr = 4.8. 

We determine the impact parameter of the encounter, 
6, randomly with probability proportional to 6, with a 
maximum, b maxi of 2000 AU. We randomize the relative 
orientations of the incoming star, the Solar ecliptic, and 
the binary orbital plane. Finally, for each incoming 
system, we run the simulation for a series of relative 
velocities between the system and the Sun. We use 
velocities of v=0. 1-2.5 km/s at intervals of 0.1 km/s 
and v=3. 0-20.0 km/s at intervals of 0.5km/s. Over 2.1 
million runs were preformed in total. 
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Fig. 2. — The effective impact parameter, b = \J (a) v /7v for exciting a Jovian to eccentricity greater than 0.1 vs. velocity of incoming star. 
They gray band indicates, given then number of trials at a given velocity, the smallest value of b we can measure with 99% confidence. 
The increase in the minimum b with velocity is a result of ours having somewhat more simulations at low than at high velocity. The right 
panel shows a zoomed in portion of the full graph shown on the left. 
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3. RESULTS 

|3.3[ we first examine the effects of 
Jovian planets, on the Kuiper Belt, 



and on Tyche in the more likely scenario where the Sun 
was bo rn in a cluster with a lifetime t\^ e = t cross . In 
Section 3.4 we examine how these results change in the 
scenario where the Sun's parent cluster was one of the 
~ 10% that reach ages of tnf e = 30 Myr. 

3.1. Effects of Close Encounters on the Jovian Planets 

Figure [I] show the probabilities of each Jovian planet 
being excited to an eccentricity greater than 0.1 as a 
function of M c and E c in a fully relaxed embedded clus- 
ter (i.e. a vir =5/3 corresponding to cr v = yj (GM C )/R C ). 
We use a range of 10 _1 -10 0,5 g cm -2 for E c because this 
covers the ran ge of e mbedded cluster properties complied 
bylFall et all (|2010|). Similarly, we use masses from 10 - 
1O 6 M because this covers the range of cluste r masses 



seen in nearby galaxies and in the Milky Way ( Lada 
Lada||2003| |Fall et al.|[2Q09| |Chandar et aL|2010p . Look- 
mg at the four graphs, we see that the probability is very 
low for all four planets in almost all combinations of E c 
and M c within the plausible range of mass and surface 
density. The probability is less than 1% for all four Jo- 
vian planets when the mass of the cluster is greater than 
1O 4 M . Neptune, the furthest planet from the Sun and 
thus the most likely to be excited, has excitation proba- 
bilities <1% over most of parameter space, and reaches 
a high of 8.7% at M c = 10 2 M©, E c = 10 5 g cm" 2 . 

A curious feature of Fig. p] is that as the mass of the 
cluster increases, the probability of an event decreases. 
To understand this, we first must realize that the num- 
ber of encounters is independent of M c at fixed E c . This 
is because A ex Uct ^ oss M c ex M°. Next, looking at equa- 
tion $h we see that as M c increases, so does the velocity 



dispersion. Figure [2] shows the velocity dependence of 
the cross section obtained though our simulation. We 
see that, at low velocities around 1 km/s, the cr oss sec- 
tion is high, compa rable to the values obtained by | Adams 
& Laughlin ( |2001 ). However, at slightly higher velocities 
the cross section dramatically decreases, dropping under 
(200 AU) 2 at around 3 km/s. Thus, at higher M c the 
velocity dispersion increases yet the number of encoun- 
ters does not. This explains the result in Fig. [I] higher 
mass clusters are less likely to increase a Jovian planet's 
eccentricity because they produce no more encounters 
(at fixed E c ), but reduce the effective cross section per 
encounter. 

Recent observational and theoretical evidence indicates 



that young embedded clusters can be sub-virial (Fiiresz 
et al.|2008NTobin et al.|2009||Offner et al.|2009| ) and thus 



we make the same calculations lor a cluster with a virial 
parameter that is one fifth the value of a fully relaxed 
cluster, a V i r =1/3. Figure [3] shows probability of excit- 
ing a Jovian planet to an eccentricity greater than 0.1 in 
a sub-virial cluster. We see that the probability is higher 
than that of a relaxed cluster, which makes sense due to 
the lower velocity dispersions, but the overall probabil- 
ity is still small. Again, Neptune is the likeliest to be 
excited with ~15% probability at the high E c , low M c 
extreme, but as in a relaxed cluster, a majority of pa- 
rameter space produces probability below 1% for all of 
the Jo vian planets. This agrees with |Proszkow &; Adams 
( 2009 ) who show that the interaction rate in a sub-virial 
cluster is greater than that in a virialized one, but not 
by much. 

3.2. Effects of Close Encounters on the Kuiper Belt 

We next check if close encounters could truncate the 
Kuiper Belt or excite an object to a Sedna-like orbit. 
To determine if the a close encounter could truncate the 
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Fig. 3. — The log of the probability of a close encounter exciting each Jovian planet's eccentricity to greater than 0.1, in one crossing 
time, as a function of cluster mass M c and surface density E c in a sub- viral cluster (a v i r =1/3). 



Kuiper Belt, we compute the expected fraction of the 
KBOs in each of our eight initial distance bins (see Sec- 
tion 2.1 for derivation). The resulting expected fractions 
in the 35-80AU and 455-500AU bins, which represent 
the extremes of our initial Kuiper Belt, are show in Fig. 
[4] for relaxed clusters and Fig. [5] for sub-virial clusters. 
We find that the expected fraction of KBOs stripped in 
the 35-80 AU bin is less than 1% for most of parame- 
ter space, only going above 1% in the high E c , low M c 
extreme, where the maximum is ~4%. The expected 
fraction of KBOs stripped in the 455-500 AU bin reaches 
88% in the high E c , low M c extreme. As we move away 
from this extreme, the probability decreases drastically. 
This appears to suggest that a birth cluster with E c and 
low M c could explain the truncation of the Kuiper Belt. 
However, comparing Fig. [4] with Fig. [I] we note that any 
cluster capable of destroying the outer Kuiper Belt while 
leaving the inner Kuiper Belt undisturbed would also 
disrupt Neptune's orbit. There is no part of parameter 
space where the Kuiper Belt is truncated but Neptune is 
not also perturbed. This is most likely because no single 
encounter can truncate the Kuiper Belt. Instead, multi- 



ple close stellar passes would be required, and the num- 
ber required is sufficiently large as to make it extremely 
likely that one of those encounters would be close enough 
to perturb Neptune. 

To see if a close encounter could produce a Sedna-like 
object, we compute the expected fraction of the KBOs 
excited to an eccentricity between 0.5 and 1.0 in a given 
distance bin. We show the result in Fig. [6] for relaxed 
clusters and Fig. [7] for sub-virial clusters. Again, we see 
that the greatest probability is found in the high E c , low 
M c extreme, with a significant drop-off as we move away 
from this part of the parameter space. Since Sedna has 
a semi-major axis of 542 AU and there is only one such 
object known, we focus on the 455-500 AU distance bin, 
and we only need a relatively small probability. We see 
reasonable probabilities exist over much of our parame- 
ter space. We conclude that it is plausible that a close 
encounter could produce a Sedna-like object. 

3.3. Effects of Close Encounters on Tyche 

To determine if Tyche could survive the birth cluster, 
we compute the cross section for Tyche becoming un- 
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Fig. 4. — The log of the expected fraction of KBOs that will 
become unbound in one cluster crossing time in a given distance 
bin. In all cases, the inner bin and the outer bin represent the 
extreme values. 
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Fig. 5. — The log of the expected fraction of KBOs that will 
become unbound in a sub- viral cluster (a v i r =1/3) in one cluster 
crossing time in a given distance bin. In all cases, the inner bin 
and the outer bin where the extreme values. 



bound. We first check to see if there was any difference 
in the cross-section depending upon which Tyche mass 
was used (see section 2.2). Using a chi-square test, we 
find that our measurements of the number of Tyche losses 
as a function of incoming star velocity at our three differ- 
ent Tyche masses are consistent with each other to ~98% 
confidence, and thus there is no reason to reject the hy- 
pothesis that the probability of Tyche becoming unbound 
is independent of it mass. Thus, we use a combination 
of data from all three masses. Looking at Fig. 
before, the highest probability, ~44%, is in the higr" 
low M c extreme, with the probability decreasing as we 
move away from this extreme. The probability is ~10% 
for moderate S c and M c . Given this, there is a small 
likelihood that Tyche would be stripped away. However, 
recall that we place Tyche at 1,000AU, and while it is 
theorized to have a semi-major axis of ~5,000 AU. Thus 
the probability of stripping likely exceeds our estimate by 
a considerable factor. Even with this increased, though, 
stripping seems unlikely in a cluster with mass > 1O 4 M . 

3.4. Encounter Effects in Long-Lived Clusters 

Finally, we consider the possibility that the natal clus- 
ter did not disperse at the point where its gas was ex- 
pelled, and instead lived on as an open cluster. As dis- 
cussed above, about 10% of embedded clusters experi- 
ence this fate. In this scenario, the Sun would be exposed 
to many more close encounters. To study the effects of 
this scenario, we use tuf e = 30 Myr in equation [3] in- 



stead of using equation [9] and repeat our analysis from 
the previous sections. 

Figure 10 shows the probability of the Jovian plan- 
ets being excited to an eccentricity greater than 0.1 in 
such a cluster. As we found for the more typical short- 
lived cluster population, the probability of disruption is 
smallest for clusters of high mass and low surface den- 
sity. Not surprisingly, excitation of the Jovians is much 
more likely in a long-lived cluster, and the only part of 
parameter space that has low probabilities is the high 
mass, low surface density extreme. Outside this regime 
all the Jovian planets are likely excited to an eccentricity 
greater than 0.1. Thus we conclude that the low eccen- 
tricities of the planets are consistent with the long-lived 
cluster scenario, but only if the cluster in question was 
both quite massive and of low surface density. 

Figures [XT] and [l2| show the expected fraction of KBOs 
lost and excited, respectively, in a long-lived cluster. We 
see that the expected fraction lost is very high for KBOs 
in all distance bins except in clusters of very high mass 
and low surface density. There is a fairly narrow strip 
of parameter space that allows the outer Kuiper Belt to 
be destroyed but leaves the inner Kuiper Belt relatively 
untouched. However, comparing to Figure [ToJ we see 
that in this parameter regime it is also very likely that 
Neptune would be driven to high eccentricity. The un- 
derlying reason is that same as for short-lived clusters: 
it takes many encounters to truncate the outer Kuiper 
Belt, and any cluster capable of producing that many 
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Fig. 6. — The log of the expected fraction of KBOs that will be 
excited to an eccentricity greater than 0.5 but less that 1.0 in one 
cluster crossing time in a given distance bin. In all cases, the inner 
bin and the outer bin represent the extreme values. 

encounters is also likely to supply one close enough to 
disturb Neptune. The behavior of Sedna is also simi- 
lar to that seen in the case of short-lived clusters: there 
is a reasonably large part of parameter space where the 
Jovian planets would not be disturbed, but there is a rea- 
sonable change of exciting a distant KBO to Sedna-like 
eccentricities. 

4. DISCUSSION 

Our conclusion that encounters in the Sun's natal clus- 
ter are likely to have very little effect on the architecture 
of the Solar System, even for very massi ve clusters, is 
at odds with much previous work, which |Adams| ( |2Q1Q[ ) 
summarizes to conclude that dynamical arguments im- 
ply that the Sun's birth cluster could not have contained 
more than several thousand stars. This conclusion, how- 
ever, is based on work that assumed the cluster lifetime 
was much greater than current observational evidence 
shows, that the surface density of a cluster varies linearly 
with the cluster's mass, and that only simulated encoun- 
ters at fairly low velocities. Each of these assumptions 
increases the probability of disrup tion encounters. For 
example |Laughlin fc Adams (19981) estimated (<r) = (230 
AU) 2 for disruption of Jovian planets due to interactions 
with binary stars in the birth cluster through scatter- 
ing simulations similar to ours. To obtain this value, 
they randomly chose the velocity of the incoming star 
from a normal distribution with a = 1 km s _1 . From 
this, they determined the rate of disruptive encounters is 
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Fig. 7. — The log of the expected fraction of KBOs that will 
be excited to an eccentricity greater than 0.5 but less that 1.0 in 
a sub-viral cluster (a v i r =1/3) in one cluster crossing time in a 
given distance bin. In all cases, the inner bin and the outer bin 
where the extreme values. 
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Fig. 8. — The log of the probability that, in one crossing time, 
Tyche becomes unbound. 



-0.13 per 100 Myr. | Adams fc Laughlinl ( [2001] ) used the 
same velocity distribution and obtained similar results. 
We can immediately identify the effects of changing ob- 
servational data which lead us to different conclusions. 
First, consulting equation ([8|, we note that a cluster of 
mass 1O 4 M and surface density 1 g cm -2 will have a 
velocity dispersion of 7 km s _1 (if it is virialized) or 3 
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Fig. 9. — The log of the probability that, in one crossing time, 
Tyche becomes unbound in a sub- viral cluster (ot v i r =1/3). 

km s _1 (for a vir = 1/5). As a result most encounters 
will be at significantly higher velocities than Laughlin & 
Adams assumed; consulting Fig. |2j we see that this will 
reduce the cross section by roughly an order of magni- 
tude. The earlier work was also based on a lifetime of 
100 Myr, while modern observations tells us that only 
~1% of stars remain in a cluste r for that long, and only 



10% survive to even 10 Myr (|Lada fc Lada 2003^ Fall 



|et al.||2009[ |Chandar et~aT7| |2Q 10] ^ X gam, switching to 
modern values leads to a 1-2 order of magnitude drop 
in the disruption probability. Together the decrease in 
cross section and in exposure time are sufficient to ex- 
plain why we find no dynamical limit on the size of the 
Sun's parent cluster. 

Changing the assumed cluster velocity dispersion and 
lifetime also explains why our conclusion s on the Kuiper 
Belt d iffers from that of earlier research. ILestrade et al.l 
( |2011| ) concluded that it was possible to truncate a debris 
disk around a main sequence star in an open cluster by a 
close encounter. To reach this conclusion, however, they 
assumed parabolic, coplanar encounters over a 100 Myr 
cluster lifetime, compared to our observationally-favored 
scenarios of mostly hyperbolic, non-coplanar e ncounters 
over a ~1 Myr c luster dissolution timescale. [Jimenez" 



Torres et al. POTT] ) created a simulation that was able to 
produce Sedna-like objects but assumed incoming star 
velocities of 1, 2 and 3 km s -1 , closest approach dis- 
tances of 200 AU and cluster lifetimes of 100 Myr. We 
concur with these authors that it is possible to truncate 
the Kuiper Belt in a cluster that lives 30 Myr, and that 
it may be possible to produce Sedna even in a shorter 
lived cluster. That the Sun might have been born in a 
long-lived cluster is by no means impossible, since ~ 10% 
of stars are. We simply point out that such a scenario 
is not typical. Moreover, we find that, regardless of the 
cluster lifetime, it is very difficult to truncate the outer 
Kuiper Belt without simultaneously disrupting the orbit 
of Neptune, a problem not considered in earlier work. 

Some other authors who have used the current obser- 
vational data about cluster lifetime and velocities have 



partially an ticipated the results presented here. Bonnell 



et al. ( 2001| ) argue that planet formation is unaffected in 
open clusters, and that it was possible for a the Sun to 



be born within a cluster whose density is greater than 
10 3 stars pc -3 if the cluster dissolves within ~10 Myr. 
Using a cluster m embership size of 1 000 and a cluster 
lifetime of 10 Myr, Adams et al. (2006) showed that typ- 
ical stellar passes are not close enough to appreciably 
enhance the eccentricity of Neptune. Their main conclu- 
sion was that planet-forming disks and newly formed so- 
lar systems generally survive thei r aggregates with little 
disruption. Spurzem et al. (2009) constructed a N-body 
simulation of the Orion Nebula Cluster, an environment 
containing massive stars similar to the hypothetical pro- 
genitors of the supernova(e) that produce the isotopes we 
see in meteorites. Their conclusion was that the critical 
threshold for the survival of wider orbit planets (those 
similar to our Jovian planets) is a cluster lieftime less 
than 10 8 yr. 



We note that |Brasser et al.| ( |2006 ) have shown through 
N-body simulations that both close encounters and clus- 
ter tidal effects have a reasonable probability of creating 
a Sedna-like object if the object begins at a semi-major 
axis >300AU, consistent with our conclusions (see Fig. [6] 
and Fig. [7]). They hypothesize a two-step formation pro- 
cess for Sedna in which the outward migration of Neptune 
scatters a KBO into an orbit with large semi-major axis, 
and the object is scattered again by passing stars. Our 
work is not inconsistent with this scenario. 

5. CONCLUSION 

In this paper we revisit the tension between meteoritic 
and dynamical constraints on the size of the Solar Sys- 
tem's parent star cluster. The meteorite evidence sug- 
gests that a supernova deposited short-lived radioactive 
isotopes near the forming Solar System. However, based 
on dynamical arguments about the likelihood of a close 
encounter disrupting the outer Solar System, |Adams fe| 
|Laughlin| ( |2001| ) and |Adams| ( |2010[ ) argue that the birth 
cluster could not have had no more than several thou- 
sand members. The probability of having a supernova 
within the range needed in a cluster of that size is only 
about 2%. Given that we have only a single example (the 
Sun) of a star system that simultaneously has supernova 
enrichment and outer planets in circular orbits, this isn't 
necessarily impossible. It could simply be that our So- 
lar System is quite unusual. However, we find that no 
such conclusion is warranted. We show that repeating 
the calculations of |Adams fc Laughlin] ( |2001[ ) but with 
embedded cluster properties drawn from recent obser- 
vations, we can update the upper limit on birth cluster 
size to values where supernova pollution is much more 
probable. 

Indeed, we show that as the mass of the birth clus- 
ter increases, the relative velocities of the constituents 
increases and the effective impact parameters for Solar 
System-disrupting events decreases, yet the number of 
expected encounters is constant. Our result differs from 
previous work because we are using current observational 
evidence that suggests that the cluster's surface den- 
sity is independent of its mass, and its lifetime is much 
shorter than previously believed. The removal of the up- 
per bound from dynamics allows the cluster size to be 
large enough to be able to produce a massive star which 
can go supernova to seed the Solar System with the ob- 
served short lived radioactive isotopes. Indeed, at cluster 
sizes of 10 4 -10 5 stars, multiple supernovae are likely, al- 
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Fig. 10. — The log of the probability of a close encounter exciting each Jovian planet's eccentricity to greater than 0.1 as a function of 
cluster mass M c and surface density S c with cluster lifetime of 30 Myr. 



lowing each supernova to be further away and exposing 
the Solar System to less mass loss. 

A close encounter within the birth cluster has also been 
proposed to explain the existence of Sedna. We find that 
this is within reason, showing that it is possible to pro- 
duce a small fraction of KBOs at ~500 AU with eccen- 
tricity between 0.5 and 1. However, we find that it is 
highly unlikely that encounters in the birth cluster could 
be responsible for truncation of the Kuiper Belt. There 
is only a small region of parameter space that allows for 
the destruction of the outer Kuiper Belt while leaving 
the inner Kuiper Belt intact, and any cluster capable 
of disrupting the outer Kuiper Belt would also be very 
likely to excite Neptune to a higher eccentricity than we 
observe. 

Finally, we caution that disruption of planetary orbits 
by the gravity of passing stars may not b e the on l y lim 



toevaporate the Sun's protoplanetary disk, preventing 
formation of the outer planets. However, this conclusion 
is highly sensitive to the rate of mass loss due to photoe- 
vaporation, which is highly uncerta in. The 10 4 star limit 



is based on a loss rate taken from Adams et al. 



but Ercolano et al. 



(2004) 



iting factor when it comes to cluster size. Adams (2010) 



argues that in a cluster of more than ~ 10 4 stars, the ul- 
traviolet radiation produced by massive stars would pho 



(2009) show that this rate is uncer- 
tain at the order of magnitude level. Moreover, it is not 
entirely clear that photoevaporation inhibits planet for- 
mation; instead, by raising the dust-to-gas ratio, it may 
trigger gravitational instability (Throop & Bally 2005). 
Due to these uncertainties, the question of whether pho- 
toevaporation might provide a limit on the cluster size 
remains an open one for future research. 
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